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ABSTRACT 

We present a theoretical study of galaxy spin correlation statistics, with de- 
tailed technical derivations. We also find an expression for the spin-density cross- 
correlation, and apply that to the TuUy galaxy catalog. The observational results 
appear qualitatively consistent with the theoretical predictions, yet the error bars 
are still large. However, we expect that currently ongoing large surveys such as 
the Sloan Digital Sky survey (SDSS) will enable us to make a precision mea- 
surement of these correlation statistics in the near future. These intrinsic galaxy 
alignments are expected to dominate over the weak lensing signal in SDSS, and 
we present the detailed algorithms for the density reconstruction for this case. 

These observables are tracers of the galaxy-gravity interaction, which may 
provide us deeper insights into the galaxy formation and large scale matter dis- 
tribution as well. 

Subject headings: galaxies:statistics — large-scale structure of universe 



1. INTRODUCTION 

The origin and evolution of the galaxy angular momentum, i.e., the galaxy spin has 
been the subject of many studies in the last century. Hoyle (1949) suggested an original 
idea that the origin of the rotational galaxy motion could be ascribed to the gravitational 
coupling with the surrounding galaxies. Sciama (1955) applied Hoyle's idea to his theory for 
the formation of galaxies in a steady-state universe model. 
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It was Peebles (1969) who first quantitatively examined Hoyle's idea in the gravitational 
instability picture. He argued that the shear effect due to the primordial tidal torquing 
from the neighbor matter distribution should be mainly responsible for the acquisition of 
the angular momentum by a proto-galaxy. He pointed out that the alternative models for 
the origin of the galaxy angular momentum such as the initial vorticity model and the 
primeval turbulence model proposed by von Wiezsacker (1951) and Gamow (1952) to a 
wrong prediction of too early formation of galaxies. Assuming a spherical symmetry of a 
proto-galaxy, he analyzed quantitatively the growth rate of the magnitude of the galaxy 
angular momentum in the frame of the linear perturbation theory, and drew a conclusion 
that the galaxy angular momentum grows as proportional to the second order perturbation 
(oc t^/^ for a Q = 1 universe). 

In contrast. White (1984) showed that the proto-galaxy angular momentum grows at 
first order (oc Hor a Q = 1 universe) unless the restrictive condition of the spherical symmetry 
is imposed on proto-galactic sites, which was originally contended by Doroshkevich (1970). 
He expanded Doroshkevich's contention in detail by means of the linear perturbation theory 
described by the Zel'dovich approximation, and confirmed that the proto-galactic angular 
momentum is generated by the misalignment between the proto-galactic inertia tensor and 
the local gravitational shear tensor, and grows to first order during the linear phase. He 
confirmed his results by N-body simulations. 

Heavens & Peacock (1988) analyzed the correlation of the galaxy angular momentum 
with the local density maxima in the linear regime, and concluded that the total angular 
momentum in the linear regime is almost independent of the height of the density peaks (see 
also Hoffman 1986, 1988). Catelan & Theuns (1996) extended the Heavens- Peacock works 
and calculated the expectation value of the angular momentum assuming an ellipsoidal 
proto-galaxy centered on a peak of the Gaussian density field using White's formula. 

While all these studies concentrated on the magnitude of the angular momentum, the 
total angular momentum, or even the fraction of virial energy in rotation, is very difficult to 
observe. On the other hand, the direction of the angular momentum, i.e., the galaxy spin 
axis can be measured only from the position angle on the sky and the projected axis ratio, 
which can be implemented for very large surveys. Therefore, the galaxy spin axis could 
provide more useful statistics that can be easily tested against real observational data. 

Recently Lee & Pen (2000, hereafter LPOO) pointed out that the Ist-order linear per- 
turbation theory predicts preferential alignments of the galaxy spin axis along with the 2nd 
principal axis of the local gravitational shear tensor, and suggested a unique statistical model 
that uses the galaxy spin axis as a tool to reconstruct the initial density field. Their theory 
is based on two basic assumptions: First, the spin axis of a galaxy aligns well with that of 
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the underlying dark halo. Second, the galaxy spin aligns with the 2nd principal axis of the 
local gravitational shear tensor to a detectable degree. 

The first assumption is generally accepted as a reasonable working hypothesis for the 
spiral galaxies. A spiral galaxy is a highly fiattened disk with its plane perpendicular to the 
direction of the underlying halo spin axis in most galaxy- formation theories (e.g., Mo, Mao, 
& White 1998). There is also an observational clue to this assumption. The galaxy spin can 
be observed to much larger radii than the galaxy radius through radio emission of the gas, 
and the spin direction of a spiral galaxy has been seen to change only very modestly as one 
moves to larger radii. This suggests that the spin axis of a galaxy is well correlated with 
that of the whole halo. 

Meanwhile the second assumption should work subject to numerical testing. In fact 
several N-body experiments have already shown that the dark halos have preferred direc- 
tion in the spin orientation, and that this preferential spin alignment has likely a primordial 
origin. Dubinski (1992) for the first time found this preferential spin alignment in N-body 
simulations of dark halos. By comparing simulations in the absence and presence of a cos- 
mological tidal field, he investigated the effect of the initial shear on the spin orientation. He 
concluded that the numerically detected preferential spin alignment measured in dark halos 
resulted from the shear effect due to the linear regime tidal torque. LPOO directly calculated 
the correlation in direction between the Lagrangian and the Eulerian angular momentum of 
dark halos measured in N-body simulations, and confirmed that the linear theory prediction 
for the orientation of the halo angular momentum is quite a good approximation. 

The central concept of LPOO is that one can use this preferential spin alignment, if 
it really exists, as a linking bridge between the initial matter distribution of the universe 
and the observable unit galaxy-spin field. They have provided a mathematical algorithm 
to reconstruct the initial shear and density fields from the observable galaxy spin axes. 
Conventionally the peculiar velocity or the weak lensing shear fields have been used to 
reconstruct the total mass density field. This new method for the density reconstruction 
using the galaxy spin field is believed to be advantageous for a couple of reasons: First, 
the galaxy spin axis is relatively easier to measure observationally. Second, it is free of the 
standard galaxy biasing. Third, it allows the reconstruction of the full three dimensional 
density field. 

Very recently, numerical study of the galaxy spin or galaxy ellipticity alignment due to 
the local gravitational shear has become quite topical. The flurry of recent activities (Croft 

& Metzler 2000; Heavens, Refregier, & Heymans 2000; Catelan, Kamionkowski, & Blandford 
2000; Crittenden et al. 2000a; Crittenden et al. 2000b) is motivated partly by the statistical 
search in blank fields for weak lensing signal, for which the intrinsic galaxy alignment plays 
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a role of systematic error. The gravitational shear effect on the galaxy spin axis due to the 
initial tidal torquing is local, distinguished from the weak lensing shear which can change 
only the apparent orientation of the galaxy spin axis due to the distant intervening matter far 
from both the source and the observer. From here on, the intrinsic cosmic shear is referred 
to as the local gravitational shear. 

Heavens ct al. (2000) calculated the intrinsic ellipticity correlation of dark halos found 
in high-resolution N-body simulations. They implied that at small redshift the intrinsic 
ellipticity correlation due to the local cosmic shear effect dominates the correlation signal due 
to the weak lensing effect. Crittenden et al. (2000a) reanalyzed the results of Heavens et al. 
and demonstrated that the results of Heavens et al. from high-resolution N-body simulations 
in fact indicate stronger intrinsic spin alignments than that of LPOO from low-resolution 
simulations. Croft & Metzler (2000) also detected the intrinsic correlation of projected 
ellipticities of dark halos in high-resolution N-body simulations, and also showed that the 
correlation signal is not strongly affected by the resolution of the simulations. Actually they 
found by comparing two simulations of different resolutions that the simulations of higher 
resolution found more intrinsic correlations. 

Observationally, the question of galaxy spin alignment has received periodic atten- 
tions. The history of observational search for galaxy alignment traces back to the 19th cen- 
tury, and has been marked by checkered records (Strom & Strom 1978; Cregory et al. 1981; 
Binggeli 1982 ; Helou & Salpeter 1982; Helou 1984; Dekel 1985 ; Lambas et al. 1988; Flin 1988 ; 
Hoffman et al. 1989; Kashikawa & Okamura 1992 ; Muriel & Lambas 2000; Codlowski 1994; 
Han et al. 1995; Cabanela & Dickey 1999). For a review of the history of the field, see Djor- 
govski (1987) and Cabanela & Aldering (1998). However, past observational searches for 
galaxy ahgnment suffered from the small sample sizes (Cabanela & Dickey 1999). It was 
only very recently that positive and reliable signals of galaxy alignments have been detected 
from large galaxy samples. Pen, Lee, & Seljak (2000, hereafter PLSOO) have reported a 
tentative detection of the intrinsic spin correlation signal from the TuUy galaxy catalog. The 
observed signal turns out to be significant at the 97% confidence level with the amplitude 
of order of 1% at Ih^^ Mpc, which is consistent with the theoretical predictions made by 
PLSOO. Brown et al. (2000) also detected the intrinsic alignment in galaxy ellipticities us- 
ing the SuperCOSMOS Sky survey data. They showed that their results agree well with 
the linear theory predictions on the galaxy preferential alignment (Crittenden et al. 2000a; 
LPOO). The observed and simulated amplitudes of correlations are expected to be stronger 
than the weak lensing effects for surveys such as Sloan Digital Sky survey (SDSS), and thus 
a quantitative analysis of intrinsic alignments must be completed before one can attempt 
to measure weak lensing shears within SDSS. These positive observational results hint at a 
possible detection of the spin-density cross correlation signal which will be addressed here. 
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The theory proposed by LPOO that the hnear shear and density fields can be recon- 
structed using the detectable intrinsic galaxy spin alignments is quite speculative, based on 
many simplifying assumptions. The idea of LPOO must go through thorough observational 
and numerical testings in the future. However, the recent observational detections of in- 
trinsic galaxy alignment and the agreement of the strength of the observed signals with the 
theoretical predictions encourage us to have a prospect for the plausibility of our theory and 
its application to the real universe. If, as predicted, the intrinsic alignment signal indeed 
dominates the weak lensing signal in shallow surveys like SDSS, the extraction of intrinsic 
shear becomes more plausible than that of the weak lensing. 

In this paper, we present the galaxy spin correlation statistics with technical details. In 
§2, we review the mathematical algorithms for the density reconstruction given by LPOO in 
greater detail for the reader's thorough understanding of our previous and future works. In 
§3 we review the spin-spin correlation statistics, and provide an analogous spin-density cor- 
relation statistics. In §4 we compare the theoretical estimates given in §3 with the observed 
signals. In §5 the results are summarized and final conclusions are drawn. We relegate the 
detailed calculations and derivations to Appendices A - J. 

2. DENSITY RECONSTRUCTION 

In the standard gravitational instability picture, a proto-galaxy acquires its angular mo- 
mentum from the local gravitational shears due to the tidal interaction with the surrounding 
matter. The angular momentum of this proto-galaxy gradually evolves till the proto-galactic 
region reaches the moment of recoUapse. On recollapse, separated out from the rest of the 
universe, its angular momentum would be approximately conserved afterwards. In other 
words, the galaxy angular momentum is expected to preserves its initial dependence on the 
local shear tensor fairly well that has been acquired during the linear regime. It is worth 
mentioning that the galaxy merging or secondary infall does not break the dependence of the 
galaxy angular momentum on the initial shears since the total rotational angular momentum 
after merging or infall process is the result of the constituent orbital angular momentum of 
the galaxies combined, which depends on the initial shear tensors. Similarly, the impact 
parameter of a collision is also determined by the shear field. Thus, what changes by those 
processes is only the smoothing scale of the intrinsic shear which correlates with the galaxy 
angular momentum. 

It is true that one cannot expect the linear theory to fully describe the evolution of the 
galaxy angular momentum. Nonlinear effects such as galaxy-galaxy interaction and etc. may 
modify the galaxy angular momentum during the subsequent evolutionary stages. Never- 
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theless recent numerical simulations have found that in fact the linear theory predictions for 
the direction of galaxy spins are in fairly good agreement with numerical results (Dubinski 
1992; LPOO). Thus, we base our study of the direction of the galaxy angular momentum on 
the linear perturbation theory. 

White (1984) and Catelan & Theuns (1996) have shown that in the 1st order linear per- 
turbation theory described by the Zel'dovich approximation, the galaxy angular momentum 
in Lagrangian space is expressed as 

Li{t) = —S'^{t) — -^-^eijkTjilik, (1) 

where S{t) is the expansion factor, D{t) describes the growing mode of the density pertur- 
bations, I = (Ilk) = (/ qiQkd^q) is the inertia tensor of a proto-galactic site in Lagrangian 
space, T = (Tji) = [djdicj)) is the local shear tensor defined as the second derivative of the 
gravitational potential, (j) smoothed on a galactic scale of R. 

Rotating the frame into the principal axis of the local shear tensor, T, we can reexpress 
equation (1) in terms of the three eigenvalues, Ai, A2, A3 of T such that 

Li OC (A2 - A3)/23, 1/2 OC (Ai - A3)/3i, L3 OC (Ai - A2)/l2, (2) 

where the three eigenvalues are ordered to be Ai > A2 > A3. We note three important 
implications of equation (2). First, if a proto-galactic region is spherically symmetric (cor- 
responding to /12 = /23 = -^31 — 0), then the region gains no angular momentum at 1st 
order. Second, if the principal axis of the inertia tensor, I, is aligned perfectly with that 
of the shear tensor, T, then the off-diagonal elements of the inertia tensor is zero in the 
shear principal axis frame, no angular momentum generated at 1st order, either. Third, 
if the proto-galactic region is non spherical, and the principal axis of the inertial tensor is 
misaligned with that of the shear tensor, then one can expect the region to acquire a net 
angular momentum vector with L2 being dominant since Ai — A3 is always bigger than the 
other two differences. In other words, the direction of the proto-galactic angular momentum 
is on average preferentially aligned with the 2nd principal axis of the shear tensor. 

For the ideal situation where the principal axis of the inertia tensor is totally independent 
of that of the shear tensor (Catelan & Theuns 1996), one can expect the maximal preferential 
alignment of the galaxy angular momentum vector along with the 2nd principal axis of the 
intrinsic shear tensor since (/I3) = (7^3) = (7^2) • What has been found in LPOO numerical 
simulations is, however, far from being idealistic. The principal axis of the inertia tensor has 
turned out to be quite strongly correlated with that of the shear tensor. However, a slight 
but detectable misahgnment between the two tensors has been detected by LPOO. It means 
that in spite of the strong correlation between the inertia and shear tensors a net nonzero 
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angular momentum at first order is indeed generated to a detectable degree with its axis 
preferentially aligned with the intermediate principal axis of the local shear. 

The essence of this idea is well represented by the following simple equation (see Ap- 
pendix A): 

{LiLj\T) = —^(^ij - (iTikTkj. (3) 

Here T is a unit traceless local shear tensor (T = T/|T| where Tij = Tij — SijTr/S), L is a 
unit galaxy spin vector, and a is a spin-shear correlation parameter introduced by LPOO to 
measure the strength of the correlation between the local shear and the galaxy spin axis. If 
a — 0, {LiLj\T) = Sij/3, spins are randomly oriented without any correlation with the local 
shears. While if the inertia and shear tensors are mutually uncorrelated, and there is no 
nonlinear effects, then the value of a is calculated to be 3/5 (it was mistakenly cited as unity 
in LPOO, see Appendix A). The real value of a should be determined empirically by numerical 
simulations, since in the linear theory one cannot estimate the strength of the correlation 
between the inertia and shear tensors from 1st principles, and one expects non-linear effects 
to be important as well. LPOO suggested the formula for the estimation of a from N-body 
simulations: 

a = 2-6A^L^ (4) 

Here {Xi}f^i are the three eigenvalues of the trace-free unit shear tensor, satisfying = 1 

and Yli Aj = 0, while L is the unit angular momentum vector measured in the shear principal 
axis frame. Note that if equation (3) holds as a theoretical estimation for LiLj , then equation 
(4) becomes optimal (see Appendix J). LPOO found a ^ 0.24 in their N-body simulations. 
For the detailed description of the measurement of a from N-body simulations used by LPOO, 
see §2 in LPOO. 

The numerical result of LPOO indicates that the present galaxy spin axes are indeed 
(weakly but detectably) correlated with the intrinsic shears even though the correlation is 
not very strong. It is worth noting that a is a universal value, independent of scale. The 
spin-shear correlation parameter of a, by its definition, must be measured from the tidal 
shears smoothed on the same scale that L is defined on. 

Given the detectable preferential alignment of the galaxy spin along the 2nd principal 
axis of the intrinsic shear, it is possible to reconstruct the shear field from the observable unit 
galaxy spins. Let us say that we have m galaxies with measured unit spins, L(x^) for 7 = 
1, 2, • ■ ■ , m. Now, we would like to find the maximum likelihood value of the traceless shear 
tensor, T at each galaxy position. Using Bayes' theorem, P(T|L) = P(L|T)P(T)/P(L). 
An immediate complication arises: P(L|T) is a purely local process, independent of the 
events at any other point, while P(T) = P[T(xi), T(x2), . . .] is a joint random process 
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linking different points with one another. As we noted in LPOO, the hnear shear expectation 
(T|L) = since P(L|T) is an even function of T. Since we are using only directions of the 
spins (which are more readily observed and predicted), we cannot recover the magnitude of 
the shear field. In other words, the shear field can be reconstructed up to the ambiguity of 
a multiplicative normalization constant. Thus we can arbitrarily normalize the shear field. 
Here we use the normalization constraint of J Tjj(x)Tij(x)(ix = 1. The nontrivial quadratic 
maximum likelihood value of the shear field with this constraint is given as the solution to 
the following eigenvector equation (Appendix D): 

J iijlm{^a,y^p)Tij{Xa)d^y^a^ -^-Tlmi^p)- (5) 

A is the largest eigenvalue of the posterior correlation operator, 

fy7m(xQ,,x^) = (Tjj(xQ,)T;m,(x/3)|L). lu the asymptotic case of a <^ 1 as in LPOO simulation 
results, iijimi'^aj^id) is given (Appendix C) as 

Cijlmi'^aj^is) — J C'ij„A;(Xo — X^)C'ijn,ofe(x^ — X^)I/„(x^)I/o(x^)fi'^X^, (6) 

where C is a two-point covariance matrix of the traceless shear tensor defined as C = 
(Cijki) = (T'y(x)T'fci(x + r)) (see Appendix B). 

In Appendix D, we explain in detail using the Lagrange multiplier method that the 
eigenvector associated with the largest eigenvalue is indeed the maximum likelihood expec- 
tation value of the shear. It is worth mentioning that in practice the posterior correlation 
function is defined only accurately at each galaxy position, so the integral in equation (6) 
must be replaced by a sum over discrete galaxy positions. We have regarded small a (a <^ 1) 
as the limit of small signal to noise. We can also find a general expression for the shear 
reconstruction in Fourier space (Appendix E): 

Here P{k) is the density power spectrum. Note that equation (7) is the optimal-filtered 
version of equation (5), holding without the constraint of small a. 

Now, the expected shear field given the unit spin field can be found as the eigenvector 
of Cijim associated with the largest eigenvalue. LPOO suggested an effective power iteration 
scheme to estimate the largest eigenvector of ^ijim- One starts with an initial guess T?-, and 
defines an iteration such that 
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^r' = , 'j-',^ +fi}-\ (8) 



Sufficiently large number of iterations converges the testing vector to the solution, i.e., the 
eigenvector associated with the largest eigenvalue with a small fractional estimation error 
proportional to (Ai/Aq)'", where m is the number of iterations, Aq and Ai are the largest 
and second largest eigenvalues respectively. Appendix F gives a general proof for the power 
iteration. 

In order to find the expected shear field by solving equation (7) using the above iteration 
method, one has to know the power spectrum of the mass density field beforehand. Here we 
describe how one can actually determine the slope of the linear power spectrum in deriving 
the shear field by equation (7): Prom an observed set of N galaxies and with an initial guess 
for the mass power spectrum, we construct a posterior shear correlation function (6), which 
is a 5A^ X 5N matrix. Prom this posterior shear correlation function, one can construct a 
weighed posterior shear correlation function given in equation (7), whose largest eigenvector 
is the shear field to be reconstructed by the iteration method described above. Here the 
largest eigenvalue A is the likelihood. We can iterate this procedure itself to measure a 
self-consistent power spectrum by varying the power spectrum to maximize A. In other 
words, at the same time when one reconstructs the initial shear field, one can also measure 
the slope of the initial power spectrum by finding such power spectrum as maximize A in 
equation (7). Note, however, that one can recover the slope but not the amplitude of the 
power spectrum since the likelihood A in (7) is independent of a multiplicative constant of 
the power spectrum. 

The final step is the reconstruction of the density field, (5(x) given the traceless shears, 
T(xj) reconstructed at each galaxy position Xj. Pirst we consider a orthonormal parametriza- 
tion of the density and the five free components of the traceless shear tensor such that 



yo = 
Vi = 



Z/2 = 



_5_ 

(-3 - V3)fu + 2V3f22 + (3 - v^)f'33 
6 

(3 - V^)^!! + 2^3^22 + (-3 - ^3)733 



6 

= V2fi2, y^^V2f23, y5^V2f3i. (9) 



y in fact is a orthonormal vector-representation of the full shear, T in terms of trace and 
traceless parts. Therefore the mutual correlation between the six components of y at the 
same position is always zero. 
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Reconstructing 5(x) given T(xi) amounts to finding (j/o(x)|yi(xi), y2(xi), • • • , y5(xi)). 
Since a linear combination of tlie Gaussian variables is also Gaussian, y is a Gaussian vari- 
able, and the covariance matrix of y, say, V {Vij = {yiVj)) can be obtained by the linear 
transformation of the shear two-point correlations (eq. [B3]). We obtain the following ex- 
pression for (7/o(x)||/i(xj),y2(xi), • • • ,7/5(xj)) (see Appendix G): 

(l/o(x)||/l(Xj),|/2(Xi),---,y5(Xi)) = J^^!^^ (10) 

<^oo 

where U = V~^, and the Greek index, v goes from 1 to 5. Equation (10) allows us to 
reconstruct the density field at an arbitrary spatial position x once the traceless shear field is 
reconstructed at each galaxy position Xj. Note that the only mathematical complication that 
arises in the reconstruction algorithm is a matrix inversion. Therefore it is computationally 
tractable, involving only linear algebra. It is worth mentioning that although the density 
field is supposed to be reconstructed in Lagrangian space, the galaxy spins are measured 
in Eulerian redshift space. We can regard this displacement between the Eulerian and 
Lagrangian spaces as noise, and convolve simply the two-point density correlation function, 
^(r) with a Gaussian filter with a peculiar velocity dispersion cr^ = 150 km/s for spiral 
galaxies (see Davis, Miller, & White 1997). 

Fig.l shows the accuracy of the reconstructed density field. It plots the cross correlation 
coefficient r = {5r{k)5m{k)) / Pr{k)Pm{k) between the reconstructed density field 5^ and 
the true density field 5m- We have implemented the algorithm with the reahstic value of the 
correlation parameter of a = 0.24 and simulated two million sample galaxy spins, and find 
a good ability of the algorithm to reconstruct the density field. As was shown in PLSOO and 
Crittenden et al. (2000a, 2000b), the intrinsic galaxy alignment is a much stronger effect 
than weak lensing for shallow surveys such as SDSS. Given that it was hoped that weak 
lensing power spectra could be measured, and thus the projected mass power spectrum, we 
would clearly expect that this stronger effect of intrinsic alignments would thus allow a better 
reconstruction of the density field in the source plane. The advantage now is that the full 
three dimensional shear field can be reconstructed, not just a two dimensional projected field 
as would be the case for the weak lensing. The above algorithm described in this section 
enables us to achieve this goal. 



3. SPIN-DENSITY CROSS CORRELATION 

In order for the algorithm given in §2 to be applied to the real universe, it is indispensable 
to have a nonzero shear-spin correlation parameter, a. Unfortunately, however, it is quite 
hard to measure the value of a directly from real observational data since it requires us to 
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Fig. 1. — Accuracy of the reconstructed density field. The sohd hne is the ideal linear case, 
while the dotted line corresponds to the strength of correlations measured in simulations and 
the TuUy data. On large scales, reconstruction is always accurate. 
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know the initial shear field beforehand (see eq. [3]). An alternative simpler way to detect 
the intrinsic spin alignment is to investigate the spatial spin-spin correlation. The local 
shear effect is due to the surrounding matter distribution. But the matter in the universe 
is spatially correlated in the standard structure formation scenario. Consequently, if galaxy 
spins are indeed correlated with the linear shear tensor by equation (3) with a nonzero value 
of a, the spatial shear correlation must induce a spatial spin-spin correlation with themselves. 

Since galaxy formation is an unsolved problem, it is difficult to make an accurate quan- 
titative evaluation of the expected level of the spin-spin correlation signal. We can make at 
most approximate analytic estimates for the order of magnitude of the spin correlation and 
its qualitative behavior. PLSOO have attempted to estimate the expected strength of the 
galaxy spin-spin correlation (see eq. [1] in PLSOO) using the Ist-order perturbation theory 
and the numerical normalization amplitude of a = 0.24. Appendix H lays out the detailed 
derivation of the spin-spin correlation function presented in PLSOO. PLSOO then measured 
the spin correlation signal directly from the observed spiral galaxies of the TuUy catalog. 
The observed signal turned out to be significant at 97% confidence level, and the amplitude 
of the signal is of order of 1% at a separation of lh~^ Mpc, in agreement with the PLSOO 
theoretical estimates. 

The consistent results of the observed spin-spin correlation signal with the theoretical 
predictions motivate us to consider a correlation of galaxy spin field with the density field. 
Galaxy spin alignment with the local gravitational shear field might result in the correlation 
of galaxy spins with the directional geometry of the nearby galaxy distribution. Perhaps 
the simplest statistic to observe is the correlation between the spin axis L and the spatial 
direction f to the nearest neighbor. Therefore, we first define a simple nontrivial spin- 
direction cross correlation function analogous to the spin-spin correlation function given by 
PLSOO such that : 

u;(r)^(|L(x).f(x)r)-u;o, (11) 

where cuq is the value of a;(r) for the case of no correlation: a;o = 1/3 for the three dimensional 
case while coq — 1/2 for the two dimensional case. Note that the two vectors, L and f, are 
both defined at the same galaxy position, x. 

For a galaxy pair at x and x -|- r, let us consider the density and shear fields smoothed 
on two different top-hat scales, say R and R', where R is the top-hat galactic radius, while R' 
is the minimum top-hat radius that encloses the galaxy pair such that R' — R + r. In order 
to avoid confusion about the smoothing scale, in this section we use an explicit notation of 
Sr and respectively to represent the density and unit traceless shear fields smoothed on 
a scale of R, while Sri and T^' for the density and unit traceless shear fields smoothed on a 
scale of i?'. A simple directional vector that one could form is the gradient of the smoothed 



-13- 



density field, V6r' (x). Assuming that galaxies form on peaks of the density field, one expects 
two neighboring galaxies to sit at the ends of a ridge connecting the two galaxies. A peak 
is by definition a location where the gradient is zero. If one considers the gradient halfway 
between the two peaks, one expects it to generically be near a saddle point, where again the 
gradient is zero. And in between, the gradient would be expected to point in a direction 
perpendicular to the peak separation r. 

Instead, one would expect that the direction of the galaxy separation vector correlates 
with the major principal axis of the local gravitational shear tensor smoothed on a scale 
of the galaxy separation. If we neglect the other two principal axes (as is often the case 
in principal component analysis), we can relate fifj to T up to a considerable ambiguity. 
But, it results in a trivial spin-direction correlation: u}{r) ~ {TiktkjTij) = due to the even 
probability distribution of T. We note, however, that the principal axis of a shear tensor is 
the same as that of its square, so we shall instead relate nfj to T^T^- analogous to equation 
(3) such that 

(r.r,|T«') = i^<5,,+6T;ff,^', (12) 

where we introduce a new quantity, a direction- shear correlation parameter of 6, to measure 
the strength of the correlation between the unit galaxy separation and the major axis of the 
local shear tensor. A careful reader may have noticed the difference of the sign ahead of the 
correlation parameters between equations (3) and (12). This sign difference arises because 
the direction of each alignment with the unit shear tensor is different, f is aligned with 
the major principal axis of {TikTkj) while L is aligned with the minor principal axis. Let 
Ai,A2,A3 be the three eigenvalues of the shear tensor, T with the order of Ai > A2 > A3. 
Then the three eigenvalues of the unit traceless shear tensor T is nothing but Ai, A2, A3 with 
= -^i — Tr/3, Tr = Y^i=i ^i- Obviously the order is the same: Ai > A2 > A3. Thus, the 
principal axes of T and T coincide. But if we consider the square of the shear, (Tj^Ifej), 
the eigenvalues are given as A^, A|, A3 but with the order of A^ ~ A3 > A|. Thus, the 
intermediate axis of the shear tensor becomes the minor axis of {TikTkj)- It explains why 
(fjfj|T) is positively proportional to TikTkj (apart from the shear- independent constant) 
while {LiLj\T) is negatively proportional to TikT^j. Also note that in equation (3) T is 
smoothed on the top-hat galactic scale of R while in equation (12) T is smoothed on the 
minimum enclosing top-hat radius of R' since the galaxy separation vector can be defined 
for a galaxy pair not for one galaxy. 

This direction-shear correlation parameter of b can be also determined in N-body sim- 
ulations in principle. We suggest the following formula for estimation of b in simulations: 

b = V2Kfl (13) 

where f is the unit separation vector in the shear principal axis frame. In practice, each is 
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obtained by measuring the separation vector of each closest galaxy pair and projecting the 
separation vector into the i-th principal axis of the local shears smoothed on the mean galaxy 
separation. Again equation (13) becomes optimal if equation (12) holds as a theoretical 
estimation formula for fjfj (see Appendix J). We have found the average value of 6 = 
0.29 ±0.01 from the same N-body simulation results that LPOO used for the measurement of 
a. It is worth mentioning, however, that the galaxy distribution is known to have a correlation 
function significantly different from that of the matter, measuring the value of b requires a 
quantitative galaxy formation model beforehand. Thus, with having no quantitative galaxy 
formation model, equation (13) provides only a qualitative approximation for the magnitude 
of b. 

With the similar method that we have used for the spin-spin correlation, one can find 
an analytic estimates of a;(r) (Appendix I) such that: 

{(jR'CrRr 

where the amplitude of A depends on the correlation parameters, a and b. It has the value 
of ab/6 and 5a6/24 for the three and the two dimensional cases respectively. Here {Sr6ri) is 
the auto correlation of the density field smoothed on two different scales of R and R', and 
aR and aji' are the corresponding rms density fluctuations: 

(5i^(x)5R,(x)) = ^ J Wth{kR)Wth{kR')P{k)47rk^dk, 

4 = (4(x)) = / WUkR)P(k)Anedk, 

4' = (4'(x)) = / W,l{kR')P{k)47rk'dk, (15) 

where the top-hat window function is given as Wth{kR) = 3[sm{kR) — kRcos{kR)]/{kRy. 

To find a closed analytic form of uj{r), we can replace the top-hat filter with the Gaussian 
filter. Using a Gaussian filter of WoikR) — exp{—k^R'^/2), and a power- law power spectrum 
of P{k) = k~'^, we find 

uj(r) = -A l^'^ (16) 

Equation (16) says that for neighboring galaxies, |ci;(r)| decreases as r~^, less rapidly than 
the spin-spin correlation that decreases as (see PLSOO). Note that in Lagrangian space 
the galaxy separation of r cannot decrease below 2>R since the top-hat radius enclosing a 
galaxy pair must be at least three galactic scale radius of R. Thus we assign uj{r) a constant 
value of uj{?>R) for r < 3i?. 
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4. SIGNAL FROM THE REAL UNIVERSE 

The unique and advantageous feature of the galaxy spin statistics presented in §2 and 
§3 is that it is a readily testable theory against real observational data since it deals not with 
the magnitude but with the direction of a galaxy spin. The spin axis of a spiral galaxy can 
be easily determined by the information of the position angle (PA) and axial ratio (7^): A 
spiral galaxy is a thin disk with a circular face-on shape, and its spin vector is perpendicular 
to the plane of the disk. Therefore, the apparent axial ratio gives the magnitude of the radial 
component of a spin vector, while the position angle determines the relative magnitude of 
the tangential components of the spin vector lying in the plane of the sky. 

In order to apply observational tests to our theory, the most suitable dataset should 
be a large sample of spiral galaxies at low redshift with the information of position an- 
gle and axial ratio (7^). The low redshift condition is required since at high redshift 
the weak lensing shear effect on the apparent orientation of the spin axis is dominant 
(Jain & Seljak 1997; Wittman et al. 2000; Heavens et al. 2000). B. Tully (2000, private 
communication) generously has provided such a galaxy catalog; The Tully galaxy catalog is 
a compilation of 35674 nearby galaxy properties over the whole sky with median redshift of 
6740 km/s. Among the total 35674 Tully galaxy properties, 12122 galaxies are identified as 
spirals. 

In measuring the spin-direction correlation, we consider all the 35674 galaxies in the 
Tully catalog to calculate the direction vectors, while we used only the spiral galaxies to 
measure the spin vectors. As mentioned in PLSOO, we suspect that the shape-shape corre- 
lation of galaxies might cause a potential problem as a form of 7?.-related systematic errors. 
The 7?.-related systematic errors are involved in the measurement of the axial ratio, TZ, found 
in the Tully catalog, caused presumably by the deviation of the shape of spiral galaxies from 
a perfect ellipse, finite thickness of galaxies, and etc. For the detailed description of the 
Tully catalog and the data analysis, see §3 of PLSOO. 

An easy way to avoid any false signal from the 7?.-related systematic errors is to measure 

the two dimensional spin-direction correlation. We project the three dimensional unit spin 
vector, L and unit separation vector, r onto the plane of the sky to obtain the two dimensional 
unit spin vector, S = S/|S|, S = L — (L ■ x)x, and two dimensional separation vector, 
t = t/|t|, t = r — (r ■ x)x. Now, the two dimensional spin-direction correlation is given by 
'-^2d('") = (|S • tp) — 1/2. Note that the projection of the spin vector onto the plane of the 
sky amounts to setting 7?. = 0, so L02D{r) is free of the 7?.- related systematic errors. 

For the three dimensional spin-direction correlation, we use an effective redistribution- 
method to deal with the 7?.-related systematic errors. We first bin the separation of every 
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galaxy pair. At each bin we uniformly redistribute the given TZ of each galaxy spin in the 
range of [0, 1) and the radial component of the given f in range of (—1, 1) (the radial direction 
is along the line-of-sight at each galaxy position). We expect that the uniform redistribution 
of TZ and f of galaxy pairs belonging to each bin eliminates effectively the systematic bias 
and false signal. Now wc rcnormalize the spin and the separation vectors after the uniform 
redistribution of TZ and r at each bin, and calculate the three dimensional spin-direction 
correlation, usoir) = (|L • f p) — 1/3. 

Fig. 2 plots the resulting observed signal (filled squares) vs. the galaxy separation, 
r — cz (km/s) with error bars, and compares the observed signals with the theoretical 
estimates. Regarding the theoretical curves in Fig. 2, we convolve the Lagrangian correlation 
(dashed fine) by a Gaussian filter with (7„ = 150 km/s to obtain the Eulerian correlation (solid 
line) the observed signal is measured in Eulerian redshift space. For a detailed description 
of the convolution procedure, see also §3 in PLSOO. The error bars are obtained from the 
experiment with the 500 sets of 12122 random two dimensional unit spins. We first generate 
the 12122 random spin vectors, and calculate the spin-direction correlation. We repeat this 
process 1000 times with different sets of random spin vectors, and compute the standard 
deviation of the spin-direction correlations. The solid line is the theoretical predictions 
given by equation (14) with the normalization amplitudes of a = 0.24 and b = 0.3 for the 
case of a power-law spectrum oi P{k) = k~'^. 

Although the observed spin-direction correlation is fairly consistent with the theoretical 
estimates qualitatively, the signal is quite weak, and the error bars arc still large. We 
expect that larger surveys like SDSS will make a precision measurement of the spin-density 
correlation signal in the near future. 

5. SUMMARY AND CONCLUSIONS 

We have presented the technical formalism in which we discuss the intrinsic galaxy spin 
correlation. We have shown how the intrinsic spin correlation is related to the initial potential 
and density fields, and how the problem can be inverted to derive the power spectrum and 
density field up to a multiplicative constant from the observable orientation of galaxy spins, as 
originally claimed by LPOO. Since the intrinsic galaxy alignments are expected to dominate 
the weak lensing signal for shallow surveys such as SDSS, our algorithm for the density 
reconstruction by the intrinsic galaxy spin alignment should be more viable than the one by 
the weak lensing shear effect. 

The formalism also allows us to address the issue of the spin-direction correlation, which 
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Fig. 2. — The spin-direction correlation signal for the case of a power-law spectrum of 
P{k) — k~'^ with the correlation parameters having the value of a = 0.24 and b — 0.3. 
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we have estimated theoretically, and measured in the observational catalog. Although the 
observed signal is reasonably consistent with the theoretical estimates, the signal is quite 
weak, and the error bars are still large. We encourage future works on the spin-direction 
correlation with larger surveys, which will make a precision measurement of the spin-direction 
correlation signal. 

We are very grateful to B. Tully for his catalog. Wc also thank U. Seljak for useful 
discussions. This work has been supported by Acadcmia Sinica and partially by NSERC 
grant 72013704 and computational resources of the National Center for Supercomputing 
Applications. 
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A. SHEAR-SPIN CORRELATION 

In order to find the expectation value of tlie unit galaxy angular momentum product 
given the unit shear tensor, (Z/jLj|T), let us first consider (LjLj|T). If the local shear tensor, 
T and the inertia tensor, I are mutually uncorrelated, the ensemble average of equation (1) 
over all orientations of the inertia tensor gives 

{LiLj\T) — eiab€jcdTakTd{IkbIld)- (Al) 

where the time-dependent proportionality constant in equation (1) is set to be unity at 
present epoch. 

From the statistical isotropy of the inertia tensor, we have the inertia tensor correlation 
ihblid) = l{5kb5id + hi^bd + SkdSbi)/^ where 7 = {{If^ + /I2 + /|3)/3) is the proportionahty 
constant. Here we stress that the density reconstruction algorithm depends not on the 
magnitude but on the direction of L. Thus, the overall proportionality constants that arise 
in the middle of our derivations can be always set to unity and we set 7 = 1 hereafter. 
Of course, the renormalized angular momentum does not have the same magnitude as the 
original angular momentum. Furthermore, it does not have the dimension of the angular 
momentum anymore. But this renormalized vector does have the same direction as the 
original angular momentum, which is all that matters. Hereafter, we will use this kind of 
renormalization frequently by setting any proportionality constant to be unity whenever it 
does not affect the direction of the angular momentum regardless of the dimensionality. 



Using (Ikblid) = {hb^id + hi^bd + hd^bi)/'^, equation (Al) can be rewritten as 

|rj-i^ ^lab^jcdiTgjPI^Qfl -\- Tg^ffTiii]^ -\- 5j^j\T\ Tik-T'kj (A2) 



On can verify that (A2) does not depend on the trace of the shear, so we rewrite it in terms 
of a traceless shear tensor Tij — Tij — 6ijTr/S such that, 

(LiLjlf) = ^Sij\f\^ - f,k%. (A3) 

Note that (LjLj|T) = |Tp for this case of independent of T and I. Since we are again only 
interested in the direction of the angular momentum vector, we may rescale equation (A3) 
to have the normalization constraint of (LjLj|T) = 1, dividing each side by |Tp such that 

2 - - 
(LjLj|T) = -^^ij — TikTkj 

= ^%+Q%-T,feTfci)- (A4) 
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The first term in the RHS of equation (A4) corresponds to the stochastic sources uncorrelated 
with the initial shear field. That is, the first term represents the ensemble average of LiLj 
for the case that the direction of the angular momentum is completely random, having no 
correlation with the shear axis, due to the modification by the other stochastic sources such 
as nonlinear effect, mutual correlation between I and T, and etc. While the second term 
in the RHS of equation (A4) corresponds to the deviation of the spin direction from the 
random average value due to it tendency to align preferentially with the intermediate axis 
of the shear tensor. 

This linear theory prediction for (LiL,:|T) holds provided that there is no nonlinear 
effects and that T and I are mutually independent. In practice, however, this condition 
is not guaranteed. Here we adopt the following simple assumption: the non-linear and 
stochastic effects are uncorrelated with the linear prediction, which adds noise to the unit 
spin vector. The nonlinear effect and the mutual dependence between T and I decreases 
the relative weight of the second shear-dependence term in equation (A4), which can be 
quantified by one parameter (say, c) such that 



As the value of c decreases, the first term of the RHS in equation (A5) dominates more, mak- 
ing the direction of the spin vector more random, the linear theory predictions with perfectly 
independent I and T corresponds to c of unity, while the completely random direction of the 
angular momentum vector corresponds to c = 0. This is the generalized quadratic relation 
we propose to express the correlation between the direction of the present galaxy angular 
momentum vector and the initial gravitational shear tensor. Retaining the framework of 
the linear perturbation theory, we treat the existence of nonlinear effect and the correlation 
between the shear and the inertia tensors as stochastic sources that tend to randomize the 
direction of the spin vector, decreasing the correlation of the angular momentum with the 
initial shear tensor. 

Now, to find the expression for {LiLj\T), let us find the conditional probability density 
function, P(L|T). The conditional probability density function, P(L|T) is usually given as 
a Gaussian distribution (see Catelan & Theuns 1996) such that 




(A5) 




(A6) 



where the covariance matrix Q is defined in equation (A5). The conditional probability 
density distribution P(L|T) can be derived by integrating out P(L|T) over the magnitude 
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of L = |L| such that 

P(L|T) = / P{L\T)L^dL = '^^^ (^L^ • Q"^ • Lj , (A7) 

where P(L|T) is in fact equal to P(L|T). Here the unit covariance matrix, Qij is given by 
equation (A5). 

In the hmit of c <^ 1, equation (A7) is simpUfied into 

P(L|f ) = i- (^1 + |[i - 3fi,nMj]) ■ (A8) 

With equation (A8) and the help of little algebra, it is straightforward to calculate the 
expectation value of {LiLj\T) in the limit of c 1 such that 



{LiLj\T) = J LiLjP(t\T)dt, 



1 c\ 3 ^ ^ 

3 + 5 ) - -^cTikTkj. (A9) 
Let us define a correlation parameter, a = 3c/5. Then, we finally get the desired expression 

(L,4|t) = ^5ij - afikfkj. (AlO) 

This equation says that for the ideal case of independent I and T, the correlation parameter 
has the value of a — 3/5 (corresponding to c = 1) while for the random spins having no 
dependence on the shear tensor a — (c — O). 

For practical purposes, it is also useful to have a similar expression to equation (AlO) 
for the two dimensional unit spins. The two dimensional unit spins mean the galaxy spins 
projected onto the plane of sky, and normalized to have a unit magnitude. Let {Si,S2) — 
(cos 0, sin 0) be the two dimensional unit spins, and P{Si, S2\T)dS — P{(f)\T)d(j) be the 
conditional probability density distribution of the two dimensional unit spins. Using the flat- 
sky approximation with L3 in the line-of-sight direction, one can say Li = -y^l — L| cos 0, 



L2 = yi-L^sin^. Then, one can say P(L I T) = P(0,L3|T). Now, P(0|T) can be obtained 
by integrating out P(0, L3IT) over L3 such that 



P(0|T) = J ^P{<P,h\l)dLs 



^|l + ~ Y + /i cos^ + /2 sin^ - 2^i2cos0sin0 j |, (All) 
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where fi = TikTki for i = 1, 2, 3, and 5(12 = TikT2k. 

Using equation (All), it is straightforward to calculate {SiSj\T) such that 



/■27r 

{Sl\f) = / cos2(/)P(0|t)#, 
Jo 



1 'in 

2--8-(/i-/2), (A12) 

f27r 



(4^|f) = / sin2(/)P(0|f)#, 
./o 



^-f(/2-/i). (A13) 

f.27r 



(^i^2|t) = / cos0sin0P(0|t)#, 



-^^12. (A14) 



B. SHEAR CORRELATIONS 

Let us calculate the spatial shear correlation, Cij^i = (Tjj(x)Tfe;(x + r)). Using Tij = 
didjcf) and = V~^(5, one can write 

Ci,ki{r) = (T,,(x)Tfc;(x + r)) 

= (6',,a,,V;25(x)a,,a,,Vr-'5(x + r)). (Bl) 

Replacing the ensemble average with the spatial average by the ergodic theorem, and 
applying the integration by parts, one can show that equation (Bl) can be rewritten as 

Cijki{v) = didjdudiV:%r) (B2) 

Using the identity relation, V^T^ = dr'il/r'"^) dr"r"'^, and with the help of little 
algebra, equation (B2) can be arranged such that 

Cijki{r) = {dijSki + SikSji + SiiSjk) { y " ^ } + in^Afi) {^i"^) + -J- ~ ~y} 

H^ifkTi + 5ikfjfi + 5iifkfj + SjkTiri + SjiVifk + hiTiTj)^^ ~ y}' ^^"^^ 

where r — r/r, Jn = nr~'^ ^{r')r'^'-^^dr'^~^. The two-point covariance matrix of the trace- 
less shears C can be also obtained by Cijki = Cijki-^kiCijnn/'^-^ijCmmki/'^ + ^ij^kiCmmnn/^- 
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C. POSTERIOR CORRELATION FUNCTION 

Let us consider the traceless posterior correlation function, 
|jj7fc(xQ,x^) = (Tjj(xQ)fl„(x£;)|L(x^)) = (T;°f;^|L'^) in the asymptotic Umit of a < 1. 
Here Xc, and x^ are two fixed galaxy positions where we would like to reconstruct the shear 
field, while x^ represents any arbitrary position of the given N galaxies such that the index 
7 = 1, 2, • • • , is dummy. Thus, it is in fact the expectation value of the quadratic shears 
given the whole galaxy spin field: 

= y"r>t^ J df'' J dff^ fl]f[l^P{T",ff^,T^\V), 

/f f Pffj'^lT" T'^ T'''! 

vf^ / dt" / df'^ T;;jf;^p(t",t^,t^)^ — ^ — ^ — — ^, 
J J P(L''') 

= J VT^ J dT'^ J dff^ f^']fi^P{f'',T^,f^)P{L'^\f^), (CI) 

where Vf^ = U^^id'i^^ P(t", t^, t-^lL^) = P(t", t^, f '^)P(L'^|t", t^, tT)/P(LT) by 
Bayes's theorem. 

Here we use the approximation of P(Lt) = P{t\ ■ ■ ■ ,t^) ^ Uj=i ^(L^) = l/(47r)^. 
This approximation can be justified as follows: Let P(L^,L^) be the joint probability dis- 
tribution of the galaxy angular momentum. One can show that this joint probability distri- 
bution can be written as P(L^, L^) = 117=1 -P(L^) + (see Appendix H). Thus, in the 
limit of a <^ 1, this approximation holds at first order of a. In this asymptotic limit, we 
also have P(L^|f t^, t^) = P(L^|t^) using that the correlation of L with t at different 
points is also C(a^) which we neglect. The constant, n^i-P(L^) — ^/i^'^)^ rescaled to 
unity in equation (CI) since the exact value of the overall constant is irrelevant to the shear 
reconstruction (any positive proportionality constant can be rescaled to unity). 

Furthermore, equation (A8) says that in this asymptotic limit of a <C 1, L-dependent 
part of P(L'^|T''') is given as 

P(L^|t^) = -af2,flLlLl (C2) 

apart from a proportionality constant (c = 5a/3). Here the L-independent part of P(L'^|T'^) 
is ignored since it does not affect the shear-reconstruction, either. Inserting equation (C2) 
into equation (CI) gives 

Wx„,x^) = -ajvT^jdT'^jdTPf^ft,f:,flP{T^,T^ 
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^ -aj d^{f^fl^J:,fl)LlLl 

Here (RyTtjfii^f^kf;!^) = J2j{T^jf2k){f!lf;!^) = J2yCijnk{^a - x^)C';„ofe(x/3 - x^) by the 
Wick theorem. We ignore the other term {T^jTim) {'^nk'^ko) since this term does not depend 
on the distance of Xq — x-j, (or x^ — x-j,), having no contribution to the shear reconstruction 
through the galaxy spins. In the continuum hmit, the sum is replaced by the integration 
over x^. 



D. MAXIMUM LIKELIHOOD EXPECTATION VALUE 

In this appendix, we provide a general argument that the maximum likelihood expec- 
tation value of a Gaussian random field can be given as the eigenvector associated with the 
maximum eigenvalue of the corresponding covariance matrix. 

Let V be a Gaussian random field. Then, the probability distribution P(v) is Gaussian 
proportional to exp(— ■ A^^ ■ v/2) where A is the covariance matrix of v. Provided that 
A is positive definite, the maximum likelihood value of v must be the one that maximizes 
exp(— v^A~^v/2) or equivalently the one that minimizes (v^- A~^-v)/2. There is an obvious 
trivial solution, v = for all points, which is of course not the solution to be sought for. 

A nontrivial solution can be found by imposing a constraint. Let us choose a quadratic 
constraint of v-^ • v = 1. Then, using the Lagrange multiplier method, we can say that the 
solution, i.e., the maximum likelihood value of v under this constraint should satisfy the 
following equation: 

where A is a Lagrange multiplier. 

Solving the above equation gives 

A-^ • V = Av, • V = 1. (D2) 

Equation (D2) says that the solution to equation (Dl) is the eigenvector of A^^ with the 
associated eigenvalue of A. Thus, the eigenvector of A~^ associated with the smallest eigen- 
value minimizes (v^- A~^-v)/2, since v^- A~-'^-v = v^-AmmV = Xmin- But, the eigenvector of 
A~^ associated with the eigenvalue of A is also the eigenvector of A itself associated with the 
eigenvalue, 1/A. Therefore, the eigenvector of A~^ associated with the smallest eigenvalue, 
Xmin is in fact the eigenvector of A associated with the largest eigenvalue, 1/Xmin = ^max- 
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Hence, one can say that the maximum hkehhood expectation value of v is in fact 
the eigenvector of the positive definite covariance matrix of A associated with the largest 
eigenvalue. 



E. INVERSION THEOREM 

This Appendix is devoted fully to prove equation (7), a nontrivial mathematical theorem 
{inversion theorem), which is at the core of our density reconstruction procedure. 

The inversion theorem says the following: If a unit galaxy spin is related to a unit 
traceless intrinsic shear tensor by equation (3) with a nonzero value of a, then it is possible 
to invert the measurable unit galaxy spin field into the initial intrinsic shear field by equation 
(7). In other words, given the unit spin field, the expected intrinsic shear field is the solution 
to equation (7) as the eigenvector associated with the largest eigenvalue of the posterior 
correlation function defined in equation (6). 

In order to prove this inversion theorem, we first prove the following three lemmas. 

Lemma 1: 

■^Im = TjjTjiTj^ =^ Aifn = (El) 

In proving Lemma 1, we do not use the Einstein summation rule, so that the repeated indices 
do not mean the summation in the following proof (but the Einstein summation rule will 
be recovered after this Lemma 1). Let us first consider the off-diagonal elements, with 
I 7^ m. 

Note that the above equation is correct only in the three dimensional case where there is 
only one choice among 1,2, 3 for the dummy index i, if i I and i ^ m. So, in the final 
term of equation (E2), the index i is not dummy. Since Tr(T) = 0, we have 

3 

^ V ^rajTjl = TiiT^i -\- Ti^Tii . (-^3) 
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Using the above equations, one can say 

j 3 

Thus, we have 

3 

/ /Tl2 I rTl2 I /Tl2 I /T^2 /Tl /Tl \ rji 

— [-'- il ~^ -'- II ~^ -'- Im ~^ -'- im ~ -'^ ii-'-mm \ J-lm 

/ rTi'2 \^ rTi'2 1^ rTi'2 rTi2 rTi2 rri rji \ rji 

~ I -'iZ ~r -^Zm -^im -^ii ^ mm ^ii-'-mm j J-lm 

= ^21^, (E5) 

since = + + 2TiiTmm and |Tp = 1 With the exactly same manner, one can also 
prove for the diagonal element. An = Tii/2. 

Lemma 2: 

Ci^Ur) = (f;,(x)T;^(x + r)) ^ C.^kiik) = {k.k, - {kikm - ^) P{k), (E6) 

where P{k) — is the density power spectrum. 
By the convolution theorem, we have 



Ci,lm{k) = T,,(k)f;^(k) = (^,k, - ^-fk'^ (^kkrn - ^-fk'^ 1$'^ 





^i3 ^2 

3 








3 J 



5lm 



(E7) 



since Tij{k) = hkj^k), Tr(T) = 5, and 5(A;) = A;2$(A;). 
Lemma 3: 



^f^T,uk) = If.Ak). m 



^''"^'^^mik) - (kk-^4]i^kikm-^-^]TUk) 



p{k) ' V 3 ; v ' 3 
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= iTiAk). (E9) 

Now, we are ready to prove the inversion theorem with the help of the above three 
lemmas. Prom here on, we will regard all the proportionality constants as unity. We will use 

equation (3) as a theoretical estimation formula for LiLj, discarding the shear independent 
Sij-term in equation (3) since it does not affect the shear inversion: 

Inversion Theorem: 

L,(x)L,(x) = -r,(x)f,,(x) =^ J |^^]|^^af(ka)^k„ = f^M, (ElO) 

iabcdi^a, ^1^^ ^ ~ J <^«w(Xa " X^jCcdji " X^)Li(x^)4(x^)(ix^. (Ell) 

By the convolution theorem, 

LbcdiK, k^) = -Cabii(K)Ccdji0^p) j LiiK + - k')4-(k')dk'. (E12) 



where 



Now, by Lemma 3, we have 



CabcdO^a, k/jj ^ 

= -y"^^%^r„,(k.) lL,{K + k^-k')L,{k')dk'dk^, 

^ -J ^pf^^ TaiK) J L,{K + k')L,{k')dk'dK. (E13) 

Let Hiji^) = Li{'x)Lj{'x) = — Tj„(x)T„j(x). Then, by the convolution theorem, Hijika + 
k^) = J Li{ka + — k!)Lj{k!)dk! . So equation (E13) is written as 

- j ^fl^fa{K)H,,{K + kp)dK = - ^pf^^f J fa{K)H,,{K + k^)dK. (E14) 

Let us define Aji{yi) = Hij{-x)fii{x) = — |T|Tj„(x)T„j(x)Tji(x). Then in Fourier space one 
can say Aji{ki3) = J Tii(ka)Hij(ka + k^)dka by the convolution theorem. Let us decompose 
A into the trace- free part and trace part such that ^jz(x) = Aji{x.) + 6jiTr{A)/3. But we 
already know from Lemma 1, the trace-free part of T'j„(x)T'„j(x)rj;(x) is given as T!,7(x) 



-28- 



(apart from the proportionality constant). Therefore, Aji{x) — ■ 
Fourier space Aji(k.p) — —Tji(k.p). Thus, equation (E14) becomes 



T|r,7(x) 





■a 



( 



) 




(E15) 



since Ccdji^ji = by equation (B3). 

But, by Lemma 3, equation (E15) is equal to Tcdikp), which finally proves the inversion 
theorem. 



We will provide a general proof for the power iteration scheme in this appendix. 

Let us assume that wc have a real symmetric positive definite n x n matrix, A, and we 
seek for the eigenvector associated with the largest eigenvalue of A. Let us say, vi, V2, ■ ■ ■ , v„ 
are the the n eigenvectors of A with the associated eigenvalues of ai, 02, ■ • • , a„ respectively 
(here we assume Oi > 02 • • • > a„ > 0). If n is not too large, then we can always find the 
eigenvectors along with the associated eigenvalues by solving the eigenvector equation, Avj = 
OiVj numerically. However, in the case n is very large, finding the maximum eigenvector by 
solving the eigenvector equation could be inefficient from a practical point of view since the 
computational time to solve the eigenvector equation could be too long. The power iteration 
scheme that wc describe and prove here is a practical method to make a fast estimate of the 
eigenvector, vi associated with the largest eigenvalue, Oi fast without solving the eigenvector 
equation for the case of large n. 

Let us start with an initial arbitrary vector, u°. Wc can construct a new vector, out 
of A and such that = Au°. Now, using the eigenvectors of A as a basis, we can expand 
u° such that u° = Yl^=i h^i- So, we can write U"*^ such that = A ^^^^=1 ^i^i — Yll=i '^M^i- 
Iterating this process m times leads to a m-th vector, u"* such that u"* — Yll=i '^i^iVi- Since 
oi is the largest eigenvalue, the first component proportional to dominates. Thus, if 
we iterate this process sufficiently large times, u"* converges effectively to the eigenvector 
associated with the largest eigenvalue. After m iterations, the fractional error caused by 
approximating Vi by u"* is proportional to (01/02)"^, which goes to zero as m becomes large. 



F. 



POWER ITERATION 
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Here, the key assumption made for this power iteration to function is that A is positive 
definite, which guarantees > for alH = 1, n. However, even for the case one has a matrix 
which is not positive definite so that not all eigenvalues are positive, one can still use the 
power iteration to find the maximum eigenvector as far as the largest eigenvalue is positive. 
It can be made by inserting secondary steps between each iteration such that 

m+l/2 

Um+1/2 _ j^^m^ ^ + u"'-\ (Fl) 

with the assumption that not all eigenvalues are negative. After the first iteration, we have 
= E^=l[l + «^/(E"=l«'^')]^'^v,. If < 0, then 1 + ai/(E;=i « -fc-) < 1- Thus, this 
refined power iteration effectively suppresses the eigenvectors associated with the negative 
eigenvalues and converges u"* to Vi. 



G. DENSITY RECONSTRUCTION ALGORITHM 

Let y = {yo,yi, ■ ■ • , y^) is a Gaussian random vector of length n + 1. The conditional 
expectation value of the first component given the rest components, (yo|yi, • " " ^Z/n) can be 
calculated such that : 



(z/olz/i, •••,?/«) = j yoP{yo\yi,---,yn)dyo, 



\yu: yi: : ynj , 

Vo —p^r. "^0 ■ 



P{yo,yi,---,yn ) 
P(yi,---,yn) 



(Gi) 



Here, since y is Gaussian, y' = (yi, • • • , y„) is also Gaussian: 



P(yo, ■■■,yn)^ ^ cxp (^ ^^'\ ''^ ) , (G2) 



V(27r)"+i|V| V 2 

where V is the (n+ 1) x (n+ 1) covariance matrix for y, and V is the nxn co variance matrix 
for y' such that V^,^ = V^^ for /i, z/ = 1, ■ ■ ■ , n, and |V|, |V'| represent the determinants of 
V, V respectively. In this Appendix, the Greek indices /i, u, r run from 1 to n. 

Let U = V~^. Then, we can rewrite 

yT.v-i.y y^-U-y 
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- Uouyuyo 



2 " " 2 

^00 / , Uo^yu\^ 1 



(yo + ^^) -7^{UooU^.-UoMou)y,yu, (G4) 
\ <^oo / ^<^oo 



But, we have: 



which can be proved by 



^ {UooU^u-Uo^Uo,) = V;,-\ (G5) 



00 



- V-^ = V, (G6) 
<^oo 



since 7^ 0, i.e., S^q — 0. 



Therefore, we can express the conditional probabihty density distribution, 

Piyo\yi, - ■ ■ ,yn) such that 

Thus, equation (Gl) can be simphfied into 

(G8) 



V27r|V||V'| 

1 ?^Oy|/i/ 



Viv||v'|-i u^i' 

To calculate |V||V'|~^, let us construct a (n + 1) x (n + 1) matrix, say V" from the 
n X n maxtrix V such that 



0, V 

Then, obviously |V"| = |V'|, i.e., \V"\-^ = |V'|-^ 

Thus one can say |V||V'|-^ = |V||V"|-^ = |V||V"-^| = |V • V"-^|. But, V • V"-^ is a 
(n + 1) X {n + 1) matrix such that 

V V""-^ =1 

V,o, I 
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where I is a n x n identity matrix. Its determinant is straightforwardly calculated to be 

|V • V"-'\ = Voo - VMrVl-;\ (G9) 
Hence, we finally find equation (G8) equal to 

(yo|yi,---,yn) = -%^, (Gio) 

<^00 

since through equations (G5) and (G9), 

C/00|V||V'|-' = C/oo(^00-Vb.yorKr'), 

= 1. (Gil) 

We emphasize that this expression for {yo\yi, - ■ • ,yn) given as equation (GIO) here is 
equivalent to equation (10) in Betschinger (1987). Bertschinger's formula is written as 
{yo\y') — ^ouV^^/yu'- Let us reexpress the matrix V such that 

Voo, Voi 

^iO 1 ^ ij 

Then, one can say ^qi'/^qq = |A*|, where A represents a cofactor of V. Now the 
determinant of the cofactor can be expressed as a sum such that |B| = ^jBojlC*! where 
Bjj = and C is the cofactor of B. This is explicitly the same as Betschinger 's formula, 
when we expand each of Betschinger's matrix elements as the determinants of cofactors, and 
sum it over Voi^. 

However, there is one advantage of our formula over Bertschinger's: Our conditional 
expectation value is often computationally much cheaper, since for a translation invariant 
random field, can be computed with the fast fourier transformation method (FFT), 
while that is not generally possible for Bertschinger's formula. 



H. SPIN-SPIN CORRELATION 



Let us first consider the three dimensional spatial spin-spin correlation, 
(|L(x) ■ L(x + r)p) = {LiL'-LjL'j). Before estimating it theoretically, it is instructive to 
understand the interpretation of L in this expression. In practice, L is the measured unit spin 
vector of an observed galaxy. In theory, however, there is no way to calculate the unit spin 
vector of an observed galaxy analytically since galaxy formation is still an unsolved problem. 
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What one can do at most (and at best) theoretically is to use some analytic estimation 
formula for the expectation value of L. Here we use equation (3) as an theoretical formula 
for L based on the linear perturbation theory. 

In equation (3) we have taken the ensemble average over the inertia tensors to obtain a 
result that depends on shear tensors. Taking the ensemble average of equation (3) over the 
shear tensors as well would result in a trivial (LiLj) — Sij/3. 

{L,L[L,L'^) = {L,L,L[L'^), (HI) 
1 + a ^ ^ ^ \ f 1 + a 

r 

= l-j + a'{f,,n,flifQ. (H3) 

We took the expectation value over the inertia tensors going from (HI) to (H2) assuming 
the inertia tensors at two different positions to be independent. 

In equation (H3) it is formidable to calculate {TikTkjT-iTij) analytically since T is 
in general not a Gaussian random field while T is. We approximate (Tj/jT)jjT/;T/^) = 
{{f,,fkjflifl^)/{mT\^)) by (f,fef,.,T;^f/.)/(|t|2)2, and apply the Wick theorem such that 



( ( - ^T^^T,, \ \ -^5,^ - oTl^Tl^ ) ) , (H2) 



{fikfkjflifL) 



rp rp rpl rpl 

J-ik-l-kj-l-ii-l-ij 

|Tp|T'p 

(|T|^|T'|2) 
9 



{T,,n,){Tl{r') + {T,,Ti){T,,T') + {T,,T'){Tu,Ti) , (H4) 



4^(0) 

since {\f\Y = 4e'(0)/9by {fijfki) = i35ik5ji+3Sii5jk-25ijSki)m/^^ (Bardeen et al. 1986). 
The r-independent term in equation (H4) is straightforwardly calculated such that 



4^(0) 3' 



(H5) 



The rest r-dependent two terms in equation (H4) can be calculated using the given Cijki 
(Appendix B). We find 

/ ir-r~, r-r~,l \ I r-r~, r-r~,l \ I r-r~, r-r~,t \ I r-r~, r-r~,t \ \ '■' 



4^2(-Q^ « \ (J/ ' \-i^j-u/ J 4C^(0) 

9 ~ ~ 

= -^^^^^ikljCkjil, (H6) 
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since Ciku — Ckjij — 0. Now using equation (B3), one can show that 

CikijCkjii = g^^ + "^"^5 ~ 4J3J5 + — J| - -Js^ + -Ja^. (H7) 

For the case of a power law spectrum of ^{r) oc r", every term in equation (H7) is proportional 
to C^(r) since J3 = ^^r" oc ^(r) and J5 = oc ^(r). Thus, one can say that for the case 

of a power law spectrum, equation (H6) is proportional to ^^(r). Let us say equation (H6) 
can be written as | + Pjr^ where /3 is a proportionality constant to be determined. Note 
that it is true only for the case of a power-law spectrum. 

In order to optimize the approximation used in equation (H4), we determine the optimal 
value of the proportionality constant by considering the real value of {TikT^jTliTlj) in the 
limit of r = 0. In the asymptotic limit of r = 0, we have 

{TikTkjTlif[-) = {fikfkjfiifij) = -{fijfij) = -, (H8) 

by equation (El). But we also have 

mkfk,f;,f;^)^^ + f3p^^. (H9) 

By equating equation (H9) with r = to (H8), we find the best-approximation proportion- 
ality constant, P — 1/6. 

So far we have treated a galaxy as a point-like object. However, a real galaxy has a 
finite size with a typical one dimensional Lagrangian scale of 0.55h~^ Mpc. Dealing with 
a real galaxy with a finite size amounts to replacing ^ by (a top-hat convolved density 
correlation on a galaxy scale of R) . Finally, we find the approximate estimation of the three 
dimensional spatial spin-spin correlation through equations (H3) to (H6) such that 

3 3 \3 6^^^ 7 
« l + J^nir), (HIO) 

with the normalization of ^ij(O) = 1. 

Note that the only approximation made in the derivation of equation (HIO) is to replace 
{fff'f') with {fff'f')/{\f\Y- We have tested the vahdity of this approximation by 
Monte-Carlo simulation, and found that it holds good for all value of r. We generated a 
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\Q-b I 1 1 1 1 1 1 — I — I— I L 

1 10 



separation (tophat units) 



Fig. 3. — The numerical verification of tlie approximation used in equation (H4). Tlie solid 
line shows the correlation function of spins on a random Gaussian lattice, while the dotted 
line is the theoretical model from equation (HIO). 
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three dimensional Gaussian shear field with the density correlation function ^(r) oc r~^, and 
smoothed it with a top-hat window on a scale radius of two grids. We then generate the 
unit spins by multiplying the shear tensors with the uncorrelated inertia tensors with an 
approximate normalization. We then plot the LHS of equation (HIO) using the randomly 
generated random fields shown as the solid line in Fig. 3, while the RHS derived from the 
density correlation function is shown as the dotted line. In this model we have used a = 3/5. 
We note the excellent agreement between the model and the simulation results, showing that 
the application of Wick's theorem to the unit shears in equation (H4) is well justified. 

With the similar manner one can also find a two dimensional spin-spin correlation, 
(|S(x) ■ S(x + r)p). For the case of a power law spectrum, all quadratic statistics in the spin 
should depend on the square of the density correlation function. It is reasonable to assume 
that the two dimensional spin-spin correlation should be also expressed as 1/2 + A^'^{r) for 
powerlaw correlations. Here 1/3 is replaced by 1/2 since (|S(x) • S(x -|- r)|^) = 1/2 for the 
two dimensional random spins. 

Through equations (A12) and (A13), we have {Sf - 5||t) = -5a{fikfki - f^kfk^)/^, 
(S'i5'2|T) = —5aTikTk2/4:. The comparison of this two dimensional expression with equation 
(3) indicates that the shear-spin correlation parameter a should be replaced by 5a/4 for the 
two dimensional case. In other words, the amplitude of A = a^/6 for the three dimensional 
case should be replaced by A — (5a/4)^/6 for the two dimensional case. Thus, we have 

^ ^ ^ ^ 1 9^ 

{S^S,S:S'^)^- + -arR{r). (Hll) 



I. SPIN-DIRECTION CORRELATION 

The same technique used in Appendix H can be also applied to calculate the spin- 
direction correlation, {LiLjfifj) such that 

= {{^-^5,j-af^,fi^l^-^8., + h (II) 

= l+^^-ahifltf^f^'fS). (12) 

Using the same approximation that is used in equation (H4) and verified by Monte-Carlo 
simulations (see Fig. 3), we have 



-36- 



Note that we calculate the shear correlation at the same position but smoothed on two dif- 
ferent scales. Following the same logic explained in Appendix H, we have 

9{fif§{fi'fif)/{Aalal,) = 1/3, 

Qi{Tifi'){fgfff) + (T;ff;f )(T5T;f'))/(4a|a|,) oc (5«(x)5^'(x))V(a|4,)- The proper- 
tionality constant can be easily obtained to be 1/6 again by considering the limit of r = 
(i.e., R' = R). 



So, finally we find 

1 ^ ab {Sr>Sr) ^ 
3 6 {aR>aR 

For the two dimensional case, replacing 1/3 with 1/2, and a with 5a/4, we also find 



(U/.r,) = , + ^^^. (14) 



1 5ab {Sr^Sr) ^ 

2 24 {aR^aR^ 



{S^S,r^r,) = - + —}^^. (15) 



J. CORRELATION PARAMETERS 

In this final appendix, we provide the derivations of the optimal estimation formula (eq. 
[4] and eq. [13]) for the two correlation parameters (a and b respectively), and the involved 
error-bar formula as well. 

Multiplying Aj to each side of equation (3), and using Aj = 0, J^i^i — Z^i — 
1/2, AjAj = —1/2, and J2i<j ■^'i -^'j ~ 1/^' equation (4) is straightforwardly derived 

from equation (3). Here note that equation (3) is used as a theoretical estimation formula 
for LiLj. The error ea involved in the measurement of the average value of a is given as the 
standard deviation of a for the case of no correlation between T and L. 

For the case of no correlation, (a) = 0. So, = -sjl^cfi). Now, by equation (4) we have 
ea = VW) = ^((2 - 6 ^ XfLlY) = (Jl) 

since (L?) = 1/3, (Lf) = 1/5 for each i = 1, 2, 3 and (LfL|) = 1/15 for each i j if L is 
random. Thus, for the Nt ensemble, we finally have = ■\jA/{bNt). 

One can derive an optimal formula for b with the similar argument. But, there is one 
notable difference in deriving the optimal formula for b. For a, we have directly used equation 
(3) to derive equation (4) where each is weighed by the square of the eigenvalue, Af . It is 
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adequate since the unit galaxy spin vector, L, is expected to be aligned with the intermediate 
principal axis of the shear tensor which does not involve any degeneracy upon a sign change 
of the shear tensor (T — > — T). Whereas, the unit galaxy separation vector, r, is expected 
to align with the major principal axis of the shear tensor, and thus orthogonal to the minor 
principal axis. But the major and minor principal axes of the shear tensor are interchanged 
by a sign change of the shear tensor. 

In order to measure the real correlation between r and the major principal axis of the 
shear tensor, we first define a secondary correlation parameter, d, by weighing each by Aj 
instead of such that 

d = Xifj. (J2) 

In the following we show that in fact b = (i/\/2, and prove that the optimal estimation 
formula for b is indeed equation (13). 

Let be the major eigenvector of T^'. Since r is supposed to be aligned with e^, 
one can model the unit galaxy separation vector as a mixture of the major eigenvector with 
a random component such that 

(r.r,|t«')=7\^ef + (J3) 

where 7^ measures the strength of the correlation between f and e^. Using = (1, 0, 0) 
in the principal axis frame of T''^', one gets {rf\T^') — I + (^ilT''^') — I ~ "3"' ^^<i 
(f||T^') = I — Now, using equation (J3) for the theoretical estimation of |fiP in equation 
(J2), we have 

, = ,^Qx._iA,-iA3). (J4) 

Approximating the three eigenvalues of the traceless shear tensor as Ai = 1/ = —A3 and 
A3 = since (Xf) = (A|) 3> (A^) and AjAj = 1, and inserting these values of the three 
eigenvalues into equation (J4), we find 7^ = dy/2. 

Now, let us model the unit spin vector as a mixture of the intermediate eigenvector of 
T^, e"^, with a random component such that 

(L,4|t^) = -«M4 + (l + a^)^, (J5) 

where measures the correlation strength between L and e^. Using = (0, 1, 0) in the 
principal axis frame of t^, one gets (L?|t^) = | + ^, (Li|f ^) = | - ^, and {Ll\T^) = 
I + Inserting the above approximate values of the eigenvalues into equation (J5), we 
have CK^ = —a/2. 
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Extrapolating equation (J3) to the limit of E! — R, let us calculate {LiLj\'T^){rifj\T^') 
at this asymptotic limit of R' = R. Using equations (14), (J3) and (J5) for the theoretical 
estimation formula, and given ■ — due to the orthogonality of the eigenvectors, we 
have 

J ir^Rv ^ ir^R'v 1 q;^7^ 1 ab 
(L,L,|T^)(f,f,|T^ )^___L^_ + _. (J6) 

It shows that b — d\f^ since o?- — —\ and 7^ = -\/2d. 

The error involved in the measurement of the average of h is also given as the standard 
deviation of b for the case of no correlation between T and f . Now we have 

y id 



Thus, for the Nt ensemble, Sf, = -^4/ (15Nt). 
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